Selective coherence transfers in homonuclear dipolar coupled spin systems 
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Controlling the dynamics of a dipolar coupled spin system is critical to the development of solid- 
state spin based quantum information processors. Such control remains challenging, as every spin 
is coupled to a large number of surrounding spins. Here we demonstrate that in an ensemble of spin 
pairs it is possible to decouple the weaker interactions (weak coupling ufj) between different pairs 
and extend the coherence lifetimes within the two-spin system from 19 fis to 11.1 ms, a factor of 
572. This is achieved without decoupling the stronger interaction (strong coupling off,) between the 
two spins within a pair. An amplitude modulated RF field is applied on resonance with the Larmor 
frequency of the spins, with amplitude u)\ , and frequency of the modulation matched to the strong 
coupling. The spin pairs appear isolated from each other in the regime where the RF power satisfies 

PACS numbers: 03.67.Lx, 76.60.-k, 82.56.-b 



Nuclear spins feature prominently in most proposals 
for solid state quantum information processors. They 
have the advantage of a simple and well defined energy 
level structure and they are normally well isolated from 
other degrees of freedom. The challenge of using nuclear 
spins in solids is to obtain control over the multi-spin 
dynamics. In a dielectric solid, the dominant interac- 
tion between the spins is the magnetic dipolar coupling. 
Since the strength of the coupling between two spins is 
inversely proportional to the cube of the distance be- 
tween them, a single spin is coupled to a large number of 
surrounding spins, and not just its immediate neighbors. 
Therefore every desired gate is embedded in a complex, 
multi-body space and its dynamics has so far proven to 
be intractable. 

Controlling the evolution of a dipolar coupled spin sys- 
tem has long been an important goal in solid state NMR, 
particularly for spectroscopic studies. For example, the 
dipolar coupling has been effectively turned off using 
techniques such as spinning the sample rapidly at the 
magic angle (9 m — cos _1 (l/v / 3)) and a variety of multi- 
ple pulse techniques, which average the spatial and spin 
tensors of the coupling respectively, as well as a combi- 
nation of these [8, 9]. 

A very useful element of control would be to map the 
physical dipolar Hamiltonian of the spin system onto 
an effective interaction that has the form of only near- 
est neighbor couplings. This would significantly simplify 
the implementation of accurate two-qubit operations in 
a many-qubit solid state spin-based quantum processor 
[1-4] . This restricted evolution is also necessary to avoid 
cross-talk between adjacent solid state quantum informa- 
tion processors in ensemble quantum computation [5, 6]. 
Without such control the gate fidelities achievable within 
a given processor element will be degraded due to leak- 
age to other members of the ensemble. Nearest neigh- 
bor mapping would also allow quantum simulations of 



many-body systems such as the Ising, XY or Heisenberg 
Hamiltonians in 1, 2 or 3 dimensions. The mapping en- 
visioned here would have significant applications beyond 
quantum information processing. For example, a nearest 
neighbor interaction could allow a more accurate deter- 
mination of distances in NMR structural studies, and 
could be used to perform sequential polarization trans- 
fers, such as along the backbone of a protein. 

Here we report the first step towards the experimental 
realization of such a scheme for the special case of an 
ensemble of spin pairs, where the dipolar coupling be- 
tween the spins within a pair is significantly larger than 
the coupling between spins on neighboring pairs. We are 
able to extend the phase memory of the spin pairs by de- 
coupling the pairs from each other, without decoupling 
the interaction between spins within a pair. The control 
sequence consists of a simple amplitude modulated RF 
field, with the modulation frequency set to the desired 
dipolar coupling strength. 

Our model system is an assembly of identical spin pairs 
with (strong) dipolar coupling w£ and weaker couplings 
between spins on different pairs. In a strong external 
magnetic field aligned along z the truncated secular dipo- 
lar Hamiltonian for this system is given by 

i Q, 0=1 i^j 

where h 1 ^ = 2ct« ct^. - ct^cA - <J^a^, and wg Q/3 is the 
coupling between spin a on pair i and spin (3 on pair j. 

The goal is to introduce a modulated RF field such 
that the effective Hamiltonian is restricted to just the 
isolated spin pairs (the first term in Eq. 1). Our solu- 
tion may be understood by viewing the RF field in the 
interaction frame of this coupling. In the fully symmetric 

case the state ct 1 + er 2 evolves as (a 1 + cr 2 ) cos (j^f-t*j + 
(o^er 2 + oier 2 ) sin (^f-t\ so we chose a RF that has 



an amplitude modulation frequency u) m — 3w£/2 and 
is given (in the lab frame) by 
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where loq is the Larmor frequency of the spins. The co- 
sine amplitude modulation produces frequency sidebands 
at ojo ± 3u)f)/2. Amplitude modulated pulses have pre- 
viously been used in NMR for simultaneously irradiating 
multiple transitions in quadrupolar spin systems [7] to 
create multiple quantum coherences in the regime where 
the RF power is significantly smaller than the strength 
of the quadrupolar coupling. 

We illustrate a simple physical picture of the averaging 
process using a 3-spin case. Consider the Hamiltonian (in 
the rotating frame) of a three-spin system, in which spins 
1 and 2 are strongly coupled and spin 3 is weakly coupled 
to spins 1 and 2 (w|> » w£J), under the RF modulation 
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The time-dependent Hamiltonian in the interaction 
frame of the (1,2) pair interaction is given by 

H(t) = e -K" s D t)h 12 /i ^ H _ ^1/^ e +i(«gt)fci 2 /4 (4) 

The zeroth order average Hamiltonian [9] of this interac- 
tion frame Hamiltonian over a period r = An/Stuf) is 



The system can then be transformed into a second inter- 
action frame via 



U' = cxp (al + o>) 



Now, in this second averaging frame of (Eq. 6), the 
first term, i.e. the residual dipolar couplings to spin 3, 
averages to zero over a cycle r' = A-k/uji. Hence the 
second averaging of the couplings to spin 3 is efficient 
when ui\ » 4w^, and the effective total system dynamics 
are generated by the (1,2) dipolar coupling (the Hamil- 
tonian of the first frame transformation). This picture 
provides the motivation for our approach, but the over- 
all dynamics are more complicated than that suggested 
by the zeroth-order average shown above. If oj^ is not 
significantly stronger than wi, the higher order terms of 
the Magnus expansion become more important. In ad- 
dition, if the strength of the 1-3 and 2-3 couplings are 



|Q1>-|10> 



|0Q>+|11> 



y |n>-|oq> . 

— |01>+|1Q> — — 



+ |co s D 



-0) 

2 D 



FIG. 1: Two pairs of strongly coupled spin-^ systems with 
each pair decomposed into its singlet and triplet manifolds (in 
the rotating frame). The triplet manifolds are weakly coupled 
to each other while the singlet manifolds do not interact. 



different, additional two and three body terms appear in 
the Hamiltonian. 

Further insight may be obtained by considering the 
energy level structure of an isolated dipole-coupled pair 
of spin- 1/2 nuclei. This has four energy levels, with 
three triplet levels corresponding to a composite 1=1 
system and a non-magnetic singlet with 1 = [10]. A 
weakly coupled set of spin pairs will largely preserve 
this structure, but transitions between the singlet and 
triplet will no longer be forbidden due to the coupling 
between spins on different pairs. Let w£J represent the 
average strength of the weaker couplings. Figure 1 il- 
lustrates how the spin pairs are decoupled from each 
other under the amplitude modulated (AM) RF irradi- 
ation. If uj m = 3uf,/2, the AM irradiation simultane- 
ously drives transitions |11) — |00) <-> |10) + |01) and 
|11)'- |00) <-» 1 00) + |11) of the triplet manifold. The 
rate at which these transitions are driven depend on the 
strength of the modulation field, uj\. If uj\ 3> tof), the 
two triplet manifolds are decoupled from each other and 
the pairs are isolated from each other. However, if the 
RF power is increased further such that u>i > cjf), the 
triplet sub-space structure gets destroyed as the strong 
coupling between the spins within the pair is decoupled. 

Thus, our scheme works in the regime where uo m — 



(6) 3w£,/2 and u% > wi > 



Not surprisingly these 



are exactly the same conditions as obtained in the previ- 
ous section. Intuitively, the RF modulation allows us to 
move into an interaction frame that is moving with the 
magnetization of the dipole coupled spin pair. The ex- 
periment bears some similarities to the spin-1 decoupling 
experiments originally proposed by Pines and coworkers 
[11, 12]. In fact they suggest that their method could 
be used to decouple a heteronuclear spin from a pair of 
identical spins. However, the scheme presented here goes 
further, and permits a coherent evolution of the isolated 
spin pair while decoupling the pairs from each other. 

It is also useful to move into the interaction frame of 
the RF modulation. The zeroth-order average Hamilto- 
nian [13] of an isolated spin pair in the interaction frame 
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FIG. 2: (a) The proton spectrum of a single crystal of gypsum 
in the [010] orientation (Pake doublet). The splitting between 
the peaks corresponds to 3w£ /2-k = 44.4 KHz. Each peak is 
broadened due to weak dipolar coupling to the other water 
molecules. Therefore, 3lo%/2it ~ 16.4 KHz (b) The narrow- 
ing down of peaks under the modulation sequence with RF 
amplitude uii = 3ojf> /4 . (c) Evolution of the single spin a x 
terms under the modulation sequence. 



of the RF modulation is 
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where the average has again been performed over one 
period of the amplitude modulation [t = 2 r K/uj m = 
47r/3a>|j). Starting from the equilibrium state where the 
spins are along the external magnetic field, a collective 
7r/2 rotation of the spins places a spin pair in the initial 
state a x + a x . This state commutes with the first term of 
the interaction Hamiltonian shown above, and the effec- 
tive evolution is only due to the second term a^a^ — a}.ai. 
The set of operators, (a x +a x , al<jy+cr y (T z , o\o 2 z — OyO^) 
form a subalgebra under the commutator that is isomor- 
phic to the Cartesian subalgebra (a x , cr y , <j z ). Thus the 
strongly coupled spins oscillate between the single spin 
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FIG. 3: (a) pulse sequence used to read out the <j\ + a x 
terms. Following the modulation pulse, a 7r/2 pulse is applied 
to rotate the a x terms to a z . During the 150 /is interval (much 
shorter than Ti) all terms other than the a z decay. A n/2 
pulse is then used to monitor a z . (b) pulse sequence used to 
read out the a^a^ + a\a\ term. Following the modulation 
two back to back 7r/2 pulses act as a double quantum filter 
to suppress the single spin a x terms. A four step phase cycle 
is necessary to implement the filter. 



state a\ + a 2 x and the two spin state cr]cr^ + o\o 2 z 
first term of the Hamiltonian in Eq. 7 were absent, this 
scheme would map onto a nearest neighbor interaction, 
and as long as the initial state of the spin pairs was within 
this subspace, leakage out of the subspace would be sub- 
stantially suppressed. However, in the current scheme, 
the initial state should be both within the subspace and 
commute with the first term of Eq. 7. For an ensemble 
of spin pairs, only the collective a x state satisfies these 
conditions. 

Gypsum (CaSC>4 -2 H 2 0) was taken as a prototypi- 
cal system for a weakly interacting ensemble of identical 
spin pairs. The protons in the waters of crystallization 
comprise the strongly coupled spins. The coupling be- 
tween protons on different water molecules is significantly 
smaller than that between protons in the same molecule. 
A unit cell of gypsum has four water molecules, with two 
pairs in two inequivalent sites. When the external mag- 
netic field is applied along the [010] orientation, the dipo- 
lar splitting at the inequivalent water sites coincide and a 
Pake doublet is observed (see Fig. 2) in a one-pulse exper- 
iment [14]. In this orientation the strong dipolar coupling 
between protons in the water molecule is w£/27t = 14.8 
kHz, and the mean coupling between protons on different 
water molecules is ujf)/2ir — 5.5 kHz. 

The experiments were carried out at room temperature 
at 7.1 T ( X H 300 MHz) using a Brukcr Avance spectrom- 
eter on a 1 mm 3 single crystal of gypsum in the [010] ori- 
entation. The length of the 7t/2 pulse used was 1.67 fis. 
The experiment was repeated as the duration of the AM 
RF was varied from 100 /is to 2.9 ms with an increment 
of 5.5 fis. The signal intensity was Fourier transformed 
with respect to the length of the modulation pulse to 
yield the spectrum shown in Fig. 2(b). Fig. 2(c) shows 
the observed ^ - <r x terms plotted against the length of 
the modulation pulse [15]. A dramatic narrowing of the 
spectral line is observed in the experiment. The effective 
T 2 of the spins under the modulation is 11.1 ms which 
corresponds to a linewidth of 29 Hz. This is a factor of 
572 times smaller than the 16.6 kHz width of a single line 
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FIG. 4: The dashed line shows the o\ + a% and the solid 
line shows the a\o-^ + g\o\ terms, (i) Under the dipolar 
Hamiltonian evolution, the above terms evolve into unobserv- 
able higher order spin correlations within 100 /is. (ii) Under 
the modulation sequence, the terms oscillate 90 deg out-of- 
phase with each other for up to 360 /j,s without any signifi- 
cant attenuation in amplitude. In this case the RF amplitude 
wi = 3wg/4. 



of the Pake doublet. 

In order to demonstrate that the spin pair continues to 
undergo a coherent evolution, we performed a second se- 
ries of experiments to specifically filter out and separate 
the u\ + a"], and the o^a\ + <j\o^ terms. The two ex- 
periments are shown in Figure 3. Figure 4(a) shows the 
coherence transfer using the dipolar Hamiltonian while 
Fig. 4(b) shows the coherence transfer under the action of 
the modulation sequence. Under the dipolar coupling the 
interactions with distant spins rapidly generate higher or- 
der spin correlations, and there is a strong damping of 
the oscillation between the single spin and the two-spin 
terms. However, under the modulation sequence this os- 
cillation is seen to extend out significantly farther. Thus 



the observed line-narrowing is not a form of spin-locking 
of the single spin terms, as occurs under strong RF ir- 
radiation, but is due to the selective decoupling of the 
weaker interactions between spins on different pairs. 

In conclusion, we have demonstrated that it is possible 
to restrict the evolution of a dipolar coupled spin net- 
work to a much smaller subspace of the system Hilbert 
space. This restriction allows us to significantly extend 
the phase coherence times for selected states. The scheme 
developed works for a system consisting of an ensemble 
of spin pairs, where the coupling between spins in the 
same pair is stronger than the coupling between spins on 
different pairs. 
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